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Abstract 

A q-discrete version of the two-dimensional Toda molecule equation is proposed 
through the direct method. Its solution, Biicklund transformation and Lax pair are dis- 
cussed. The reduction to the q-discrete cylindrical Toda molecule equation is also dis- 
cussed. 
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1. Introduction 



The discovery of the quantum groups[l,2] gave rise to a new phase to the studies of 
integrable systems. It also shed hght to the theory of so-called q-analog of the special 
functions, where the q-difference operator plays a similar role to the differential operator 
in the theory of the special functions [3-5]. On the other hand, it is known that the soliton 
equations have a close relationship with the special functions[6,7]. Hence it arises as an 
interesting problem to study "q-discrete integrable systems" . 

In a previous paper [8], we have proposed the q-discrete version of the two-dimensional 
Toda lattice equation and shown that the q-discrete cylindrical Toda lattice equation is 
derived through a reduction. Moreover, we have shown that its solution is expressed by a 
determinant whose entries are the q-Bessel function. In this letter, we will discuss the q- 
discretization of the two-dimensional Toda molecule (2DTM) equation (finite, non-periodic 
two-dimensional Toda lattice) [9] and derive its solution, Backlund transformation and Lax 
pair. We also present the q-discrete cylindrical Toda molecule equation. 

Let us first give a brief review of the 2DTM equation, 

^^ur,-^-UM _^UM-uj,+^ iV = l,...,M, (la) 

oxoy 

Uq = -OO, Um+1 = +00 , (lb) 

which is rewritten as 



dV 

= Vn{Jn - Jn+i) , (2a) 

^ = Vtv-i - VW , (2b) 
dy 

Vo = Vm = 0, (2c) 

by a suitable variable transformation. Equations (2a)- (2c) are transformed to the bilinear 
form. 
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-Tn 



Otn dr 



N 



dxdy dx dy 

through the dependent variable transformations, 



— tn+iTn-i 



(3) 



Vn 



Jn = ^ log 

ox Tn 



(4) 



The solution of the bilinear form (3) is given by the two-directional Wronskian, 



f 
dyf 



dxdyf 



In (5), / is chosen to satisfy the boundary condition (2c) as 



(5) 



M 

f ^^(t>k{x)xk{y) , (6) 
fe=i 

where (pkix) and Xkiv) are arbitrary functions in x and y, respectively. 

Our method relies on the fact that the bilinear forms of "generic" soliton equations are 
nothing but the identities of the determinants or the Pfaffians, or those derived from their 
reduction. In the present case, (3) is reduced to the Jacobi identity of the two-directional 
Wronskian. This fact can be a powerful guiding principle for the extension of the integrable 
systems. For example, usual discretization of soliton equations is performed based on the 
principle [10]. We will apply the idea to the q-discretization of the soliton equations. 



2. q-Discrete 2DTM Equation 



We propose a system given by 

Sq'',xVN{x,y) = VN{q°'x,y)JN{x,q^y) - yjv(x, y) Jjv+i(x, y) , 
5q0^yJN{x,y) = VN-i{q"x,y) - VN{x,y) , 
Vo{x,y) = VM{x,y) = , 
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(7a) 
(7b) 
(7c) 



where Sqa ^. and 5qi3 y are the q-difFerence operators defined by 

The operators Sqa ^^ and (^g/s^y tend to and in the hmit — > 1, respectively. In 
this hmit, (7a)-(7c) are reduced to the 2DTM equation (2a)-(2c). We caU eqs. (7a)-(7c) 
the q-discrete 2DTM equation. 

Equations (7a)- (7c) are transformed into the bihnear form, 



V,a; SqP^yTN{x, y) ■ Tn{x, y) - Sq^^xTNix, y)5q0^yTN{x, y) 

= Tn+1 {x, y) Tn-1 (q^'x, ql^y) , (9) 



through the dependent variable transformations. 



JNix,y) 
VNix,y) 



TN-i{x,y)TN{q°'x,y) 
(1 - q)x {TN-iiq'^x,y)TNix,y) 

TN+i{x,y)TN-i{x,q^y) 
TN{x,y)TN{x,q^y) 



- 1 



(10a) 
(10b) 



The solution of the bilinear form (9) is given by the two-directional Wronski-type deter- 
minant. 



Tjv(x,2/) 



f{x,y) Sq^,^f{x,y) 

Sq0,yf{x,y) Sqc.^jq0^yf{x,y) 



Sq^,xf{x,y) 



(11) 



<5f,-V(x,y) Sq.,J^,-;f{x,y) ■■■ <«7J<5j-V(x,y) 
where f{x,y) is chosen to satisfy the boundary condition (7c) as (6). 

Let us prove that (11) really gives the solution of the bilinear form (9). We have for 
example. 
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f{q '^x,y) Sq^,^f{q '^x,y) 



f{x,y) Sg<.^a:f{x,y) 

^q'3,yf{x,y) Sqc.^Jgi3yf{x,y) ■ 



8q'-~x^qP,yf{(r'^X^y) 



CPqi',yfir^^^V) 



S^.Jf{q--x,y) 

^q-~x^qP,yf{(l~'^X^y) 



xN-1 xN-1 

0„a ™ 



fi<l '^x,y) 



(12) 



In deriving the second determinant from the first, we have subtracted {k + l)-th column 
multiplied by {1 — q)q~°'x from A;-th column for fc = l,---A^ — 1, to confine the shift of the 
independent variable x to the most right column of the determinant. Moreover, we note 
that S^cT^ f{q~°-x,y) means 5^J'^^ f{x,y)\x^q-ay;. Multiplying N-th. column by {l — q)q~°^x 
and adding (A^ — l)-th column to A^-th column in the second determinant of (12), we get 



(1 ~ q)q riv(g "x,|/) 



fix.y) 



f{x,y) 

3q^,yf{x,y) Sq°',xSqP,yf{x,y) 



q' ,y 



f{x,y) Sqc^^^S^^ ^f{x,y) 



q'^,y 



Similarly, we obtain 



(1 ^ g)2g-(«+^)xy TN{q~'^x,q-^y) 



f{x,y) 



tN-2 



f{x,y) 



^q.-:yfM 

S^,-^'f{x,q-^y) 



C'^fiq-'-x.y) 

Sqc~xSqP,yf{q~'^X,y) 



S^^-xX^:y'fix,y) 



^q",x 



s':^-:xX,-:y'fiq-x,y) 



q' ,y 



Applying Jacobi's identity on (14) with A*" replaced by A/" + 1, we obtain 



(13) 



(14) 



TN{q "a;,g ^y)TN{x,y) - TN{q °'x,y)TN{x,q ^y) 
= {l-q)\-^^^^^xy TM+i{q-''x,q-f'y)TM-i{x,y) , 
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(15) 



which is nothing but the bihnear form (9) with x and y replaced by q~°'x and q~^y, 
respectively. Thus we have proved that (11) gives the solution of (9). 

We now discuss a reduction of the q-2DTM equation. Putting xy = and a = P = 2, 
and imposing the condition that tn{x, y) depends only on r, we find that the bilinear form 
(9) and its solution (11) are reduced to 

1 2 
{-5q,r + q5l^r)^N{r) ■ TN{r) - {Sq,rrN{r)} = TN+i{r)TN-i{q'^r) , (16) 

and 



rjv(r) = q-^N{N-i){N-2) ^-n{n-i) 

rSq,rf{r) {r5q,r?f{T) 



(r(5,,,)^-V(r) 
{r5q,rrf{r) 



(17) 



respectively. Equation (16) tends to the cylindrical Toda molecule (cTM) equation[ll], 



(18) 



in the limit q 1, and hence we call (16) the q-cTM equation. Note that (16) is trans- 
formed to 



5g,rViv(r) = qJN{qr)VN{qr) - JN+i{r)VN{r) , (19a) 
{q5q,r + l)JNir) = VN{r) - VN-i{qr) , (19b) 
Voir) = FM(r) = , (19c) 



through the dependent variable transformations, 
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rN-i{qr)TN+iir) 
rN{r)TN{qr) 
1 { rM{qr)TN-i{r) 



(1 - q)r \ TN{r)TN-i{qr) 



- 1 



(20a) 
(20b) 



3. Backlund Transformation and Lax Pair 

By using the fact that the solution of the q-2DTM equation (7) is given by (11), we 
here propose the Backlund transformation. It is written by 



= -Tjv+i(a;,y)T^_i(a;,g^y) , 
Sq^,x TN{x,y) ■ TN_i{x,y) - TN{x,y) 5qc^^ r'^-iix^y) 
= TN-iiq'^x,y)T'^{x,y) , 



which transforms a solution of the q-2DTM equation. 



TN{x,y) = 
to another solution. 



f{x,y) Sq^^^f{x,y) 

Sq0,yf{x,y) Sqc^^5q0^yf{x,y) 

S^,-^'f{x,y) 5q.,J^i-2f{x,y) 



Sq-',xf{x,y) 



<«7.''^J:;/(^'2/) 



(21a) 



(21b) 



(22) 



T'Nix,y) 



dq^,xf{x,y) 5l'-,xf{x^y) 

^q°',xSq0,yfix, y) 5^o. Jq0^yf{x, y) 



5!" 



g" ,x 
'q",x^ql^,y 



f{x,y) 



Sq^,xS^,-;f{x,y) 5l.J^,-^f{x,y) ■ 
In other words, (21a) and (21b) are the identities for the determinants (22) and (23). This 
fact is shown by the Plucker relation as follows. 



(23) 



Let us prove the second equation (21b). First, we introduce notations 



tn{x, y) = \0,1, N -1\ , 



T'j^{x,y) = \l,2,---,N\ . 



Namely, the number "/c" in (24) and (25) means a column vector 



Then we have 



and 



> N 



TN-i{x,y) = |O,l,---,iV-2,0| , 



TNiq-'^x, y) = |0, 1, • • • , TV - 2, AT - 1 , 



r'N-i (9""^, 2/) = 1 1, 2, • • • , AT - 2, iV - 0| , 



(1 - q)q-''x T'j,{q-^x, y) = \1,2, ■ ■ ■ N - 1, N - 1 



where 



"AT - 1 



q~"'x 



( <«7x/(9-"a;,2/) \ 

^q'-~x^q0,yf{q~'^X,y) 

\5^.-^5^,-^f{q-x,y)J 



and 









\1/ 



(32) 



which is inserted to equahze the size of the determinant. 
We now consider an identity of 2N x 2N determinant, 







Oil ••• N-2 \ 

i i 1 







N-l N-1 



q-°'x 







N -2\ N -I N-l 



(33) 



Applying the Laplace expansion to the right-hand side, we obtain an identity ( the Pliicker 
relation) , 



0=|0,l,---,7V-2,Ar-l,-„,| |1, • • • , TV - 2, TV - 1, 0| 
-|0,l,---,Ar-2,Ar-l| |l,---,iV-2,7V- V«^,0| 
-|l,---,iV-2,iV-l,iV- V«,| |O,l,---,iV-2,0| , 



(34) 



or equivalently. 



T7v(g ''x,y)r^_i(a;,y) - rjv(a;,y)r^_i(g "'x,y) 

= (1 - q)Q~'^x Tjv-i {x, yy^iq-'^x, y) , (35) 

which is nothing but (21b) with x replaced by q~°'x. Thus we have completed the proof. 
The first equation (21a) is proved in a similar way. 

It is possible in general to construct the Lax pair from the Backlund transformation. 
Following the method developed by Hirota et. al. [12] , we derive the Lax pair for the q-2DTM 
equation (7) from (21). Introducing •0 by 

TNi^i y) = TNix, y)il)N+i{x, y) , (36) 
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we have from (21), 



Sq^,xipN{x,y) = -JN{x,y)ipN{x,y) - ipN+i{x,y) 



(37a) 
(37b) 



Let us define two matrices L and R by 



/ 

Viix,y) 



L{x,y) 



\ 



\ 







VM-i{x,y) 0/ 



(Ji{x,y) 1 \ 
J2{x,y) 1 



R{x,y) = - 







JM-iix,y) 1 

JMix,y)J 



Then (37) are rewritten as 



where 



(38) 



(39) 



Sg0^y^{x, y) = L{x, y)^{x, q^y) , Sq.,^^{x, y) = R{x, y)^{x, y) , (40) 



/ "^lix^y) \ 
^{x,y)= : 

VV'M(a;,y)/ 

The compatibihty condition of the hnear system (40) yields 



(41) 



5q.^^xL{x, y) - Sq0^yR{x, y) = R{x, y)L{x, y) - L{q'^x, y)R{x, q'^y) , (42) 

which recovers the q-2DTM equation (7). Consequently, (38) and (39) give the Lax pair 
of the q-2DTM equation. 
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4. Concluding Remarks 

In this letter, we have proposed the q-2DTM equation, its solution, Backlund trans- 
formation and Lax pair. There are several definitions of integrability for the continuous 
equations, such as the existence of N-soliton solution, that of an infinite number of con- 
served quantities or symmetries, Painleve property and so on. For a given continuous 
equation, there are several ways to discretize it, depending on the definition to be taken. 
In our case, the guiding principle of discretization is to preserve integrability in such a 
sense that the determinant structure of solutions of the discretized equation is the same 
as that of the original continuous equation. For the discretization we have employed the 
bilinear formalism, since the determinant structure is clearly seen in it. 

It was revealed that the time evolution of the solutions of the usual integrable system 
is subject to the infinite dimensional Lie algebra[13]. A question naturally arising is what 
the algebra describing the q-discrete system is. Is it the quantum group or some others? 
It may be an interesting and important problem to find the structure of the algebras 
underlying the q-2DTM equation. 

Finally, we mention the possibility of q-discretization of the Painleve equations. There 
is a close relationship between the Painleve equations and Toda molecule equation. 
Okamoto has shown that the r functions of Painleve equations satisfy the several types of 
Toda molecule equations[14]. In particular, the r function of the third Painleve equation 
satisfies the cTM equation. Moreover, it has been shown that the Painleve equations are 
reduced to the identities of determinants[15,16] for special values of the parameters in the 
equation. Since we have the q-cTM equation, it is natural to expect that we can perform 
the q-discretization of the third Painleve equation. We will report on this subject in a 
forthcoming paper. 

References 



11 



[1] V. G. Drinfel'd, Sov. Math. Dokl. 32(1985) 254. 
[2] M. Jimbo, Lett. Math. Phys. 10(1985) 63. 

[3] H. Exton, q-Hypergeometric Functions and Apphcations (EUis Horwood, Chi- 
chester, 1983). 

[4] T. H. Koornwinder, in: Orthogonal Polynomials ed. P. Nevai (Kluwer Academic, 

Dordrecht, 1990) p.257. 
[5] T. Masuda, K. Mimachi, Y. Nakagami, M. Noumi and K. Ueno, J. Func. Anal. 99 

(1991) 357. 

[6] A. Nakamura, Prog. Theor. Phys. Suppl. 94 (1988)195. 

[7] K. Okamoto, in: Algebraic Analysis (Academic, Boston, 1988) p. 647. 

[8] K. Kajiwara and J. Satsuma, J. Phys. Soc. Jpn. 60(1991) 3986. 

[9] A. N. Leznov and M. V. Saveliev, Physica 3D (1981) 62. 
[10] Y. Ohta, R. Hirota, S. Tsujimoto and T. Imai, submitted to J. Phys. Soc. Jpn. 
[11] R. Hirota and A. Nakamura, J. Phys. Soc. Jpn. 56(1987) 3055. 

[12] R. Hirota, S. Tsujimoto and T. Imai, to appear in Proc. of Future Directions of 

Nonlinear Dynamics in Physics and Biological Systems (1992). 
[13] M. Jimbo and T. Miwa, Publ. RIMS, Kyoto Univ. 19(1983) 943. 
[14] K. Okamoto, Ann. Mat. Pura Appl. 146(1987) 337; Japan J. Math. 13(1987) 47; 

Math. Ann. 275(1986) 221; Funkcial. Ekvac. 30(1987) 305. 
[15] A. Nakamura, J. Phys. Soc. Jpn. 61(1992) 3007. 
[16] K. Kajiwara, Y. Ohta and J. Satsuma, in preparation. 



12 



